A systematic study of neutron magic nuclei with N = 8, 20, 28, 50, 82, and 126 

in the relativistic mean field theory 
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We perform a systematic study of all the traditional neutron magic nuclei with N = 8, 20, 28, 50, 
82, and 126, from the neutron drip line to the proton drip line. We adopt the deformed relativistic 
mean field (RMF) theory as our framework and treat pairing correlations by a simple BCS method 
with a zero-range 5-force. Remarkable agreement with the available experimental data is obtained 
for the binding energies, the two- and one-proton separation energies, and the nuclear charge radii. 
The calculated nuclear deformations are compared with the available experimental data and the 
predictions of the FRDM mass formula and the HFBCS-1 mass formula. We discuss, in particular, 
the appearance of sub-shell magic nuclei by observing irregular behavior in the two- and one-proton 
separation energies. 



I. INTRODUCTION 

"Magic number" is a very important concept in many 
subjects of physics, such as atomic physics, nuclear 
physics and micro cluster physics. In nuclear physics, nu- 
clei with magic numbers have been hot topics of nuclear 
research since the beginning of this subject 0,0- In re- 
cent experiments with radioactive nuclear beams (RNB), 
disappearance of traditional magic numbers and appear- 
ance of new magic numbers are observed in nuclei with 
exotic isospin ratios. Furthermore, nuclei with magic 
neutron numbers in the neutron-rich region are of spe- 
cial interest for the study of astrophysical r-process Q ■ 

The unusual stability of nuclei with neutron (proton) 
numbers 2, 8, 20, 28, 50, 82, and 126, commonly re- 
ferred to as "magic numbers" , was traditionally explained 
in the nonrelativistic shell model approximated by the 
3D Harmonic-Oscillator central potential together with 
a very strong spin-orbit interaction introduced by hand 
On the other hand, in the models within 
the relativistic framework, say the relativistic mean field 
(RMF) theory pj IE IE HE LLif i the strong spin-orbit in- 
teraction appears naturally as the interplay between the 
strong scalar and vector potentials, which are necessary 
for reproducing the saturation properties of nuclear mat- 
ter. Due to the proper setting of the scalar and the vector 
potentials the shell structure is obtained without any ad- 
ditional parameters for the spin-orbit splittings. 

Recently, it is argued that the magic numbers are of 
a localized feature, i.e. in nuclei with exotic isospin ra- 
tios the classical magic numbers do not necessarily hold 
and new magic numbers may appear. There have been 
several experimental evidences supporting such a belief, 
including the lately reported N = 82 shell quenching 
|12| and the appearance of a new neutron magic num- 
ber N = 16 in the neutron- rich light nuclei [Rj- Such 
a localized feature of magic numbers are also claimed to 
be important for various nuclear-astrophysical problems 
[l4|. The disappearance or quenching of nuclear magic- 
ity near both neutron and proton drip lines have been 
discussed quite a lot in various nuclear models, includ- 



ing the infinite nuclear matter (INM) model [l5j , the ex- 
tended Bethe-Weizsacker mass formula Il6l , the antisym- 
metrized molecular dynamics (AMD) [13, the Hartree- 
Fock method p~3] > an d the relativistic mean field (RMF) 
[l^.l20l | theory. In the relativistic mean field model, while 
the Z = 8, 20, 28, 50, 82, and 126 isotopic chains have 
been discussed a lot, a systematic study of all the nuclei 
with neutron magic numbers is still missing. 

The appearance and disappearance of the magic num- 
ber effect are intimately related with the spin-orbit in- 
teraction and the onset of deformation. In this sense, 
a systematic study of nuclei in various mass regions in 
terms of the relativistic mean field theory is suited, since 
the spin-orbit interaction arises naturally in relation with 
the saturation property. Hence, we are able to provide 
the variation of the spin-orbit splittings with the pro- 
ton and neutron numbers. As for the deformation, nuclei 
with magic numbers are not easy to acquire deformation, 
because the onset of deformation does not utilize the en- 
ergy gain coming from the large shell gaps. At the same 
time, due to this reason it is easier to find the appearance 
of the sub-magic effects for nuclei with either the proton 
number or the neutron number being magic numbers. 

In the present work, we would like to make a systematic 
study of neutron magic nuclei from the neutron drip line 
to the proton drip line in terms of the deformed relativis- 
tic mean field theory. The pairing correlations are also 
very important to make the nucleus tend to be spherical. 
Therefore, we take the recently proposed method of using 
the zero-range (5-force for the pairing correlations in or- 
der to pick up the resonant states in the continuum when 
the nucleus approaches the neutron threshold 21]. The 
quadrupole constrained RMF calculation is performed for 
each nucleus to find all the energy minima as a function 
of the nuclear deformation. In this paper, we stay our 
calculations at the mean-field level and hence anticipate 
that we may miss the onset of small deformation in this 
framework for the case, in which the deformation effect 
is not dominant. Hence, the purpose of this paper is 
to make a systematic study of the neutron magic nuclei 
in the mean field framework and compare with the avail- 
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able experimental data in the global sense and with other 
more sophisticated theoretical models for particular nu- 
clei. We mention that even at the mean-field level, many 
proton magic nuclei appear to be deformed |20j. 

This paper is organized as follows. We provide a short 
summary of the deformed RMF+BCS method in Sec. II. 
In Sec. Ill- VI, we present the results of our calculations, 
including the binding energies, the two- and one-proton 
separation energies, the nuclear radii, and the nuclear 
deformations. With these quantities, we also discuss the 
magicity of these neutron magic nuclei with N = 8, 20, 
28, 50, 82, and 126. The whole work is summarized in 
Sec. VII. 

II. THEORETICAL FRAMEWORK 

In the present work, the recently proposed deformed 
RMF+BCS method with a zero-range (5-force in the pair- 
ing channel is adopted [2l|. The zero-range (5-force has 
proved to be very successful to take into account the con- 
tinuum effect by picking up the resonant states both in 
relativistic and nonrelativistic self-consistent mean field 
models [lH El |H |H IH El ■ In the mean-field part, 
the TMA parameter set is used j27|. The RMF calcula- 
tions have been carried out using the model Lagrangian 
density with nonlinear terms for both a and uj mesons as 
described in detail in Refs. 0,123, which is given by 

- -Sl^W + \ m l^^ + ig 4 (+vc^) 2 - g^^iptJn 

- iF^-efoJ-p-A^, (1) 

where all symbols have their usual meanings. The corre- 
sponding Dirac equations for nucleons and Klein-Gordon 
equations for mesons obtained with the mean-field ap- 
proximation are solved by the expansion method on the 
widely used axially deformed Harmonic-Oscillator basis 
[2lL l28| . The number of shells used for expansion is cho- 
sen as Nf = Nb = 20. More shells have been tested for 
convergence considerations. The quadrupole constrained 
calculations have been performed for all the nuclei consid- 
ered here in order to obtain their potential energy sur- 
faces (PESs) and determine the corresponding ground- 
state deformations |2ll |29| | . For nuclei with odd number 
of nucleons, a simple blocking method without break- 
ing the time-reversal symmetry is adopted 0, l3(l | . The 
pairing strength Vo is taken to be the same for both pro- 
tons and neutrons, but optimized for different regions by 
fitting the experimental two- and one-proton separation 
energies, more specifically, for the N = 8, 20, 28, 50 iso- 
tonic chains, V = 344.1 MeV fm 3 ; for the N = 82, 126 
isotonic chains, Vq = 310 MeV fm 3 . 



Whenever the zero-range 6 force is used either in the 
BCS or the Bogoliubov framework, a cutoff procedure 
must be applied, i.e. the space of the single-particle 
states where the pairing interaction is active must be 
truncated. This is not only to simplify the numerical 
calculation but also to simulate the finite-range (more 
precisely, long-range) nature of the pairing interaction in 
a phenomenological way |3ll |32( . In the present work, 
the single-particle states subject to the pairing interac- 
tion are confined to the region satisfying 

e l -X<E cut , (2) 

where Ci is the single-particle energy, A the Fermi energy, 
and E cut = 8.0 MeV. We find that increasing _E cut from 
8.0 MeV up to 16.0 MeV, followed by a readjustment 
of the pairing strength Vq, does not change the results 
appreciably, and therefore none of our conclusions will 
change. Recently, such a cutoff issue has been discussed 
in much detail by Goriely et al. in the Hartree-Fock- 
Bogoliubov (HFB) framework [3^. Our cutoff procedure 
happens to coincide with their conclusions (see Table II 
of Ref. although a bit different from their final 

optimal choice used in the construction of the HFB-2 
mass table. 

Even after the pairing interaction between like nucle- 
ons is taken into account, mean-field calculations (HF or 
RMF) still miss some residual correlations, one of which 
is the so-called Wigner effect. Such an effect is known 
to be important for light N * Z nuclei by making them 
more bound by about 2 MeV. Although it is usually ap- 
proximated in a phe nomenological way in most mass ta- 
bles [H |H |H HI, here we will ignore this effect 
due to the following considerations. It is expected that 
these residual correlations (including the Wigner effect) 
have only small contributions and are active only under 
special situations. In the case of the Wigner effect, it 
becomes appreciable only for light nuclei with N « Z. 

Finally, the center-of-mass correction is approximated 

by 

E cm = -J41A- 1 / 3 , (3) 

which is often used in the relativistic mean field theory 
among the many recipes for the center-of-mass correction 
|37| . Since adopting new schemes implies that one has 
to readjust the model parameters, we will adhere to this 
simple approximation in the present work. 

III. BINDING ENERGIES 

The nuclear binding energy is one of the most basic 
properties of nuclei. Although the two- and one-proton 
separation energies are more useful for the purpose of 
identifying the shell structure, the binding energies of 
nuclei can also show some important features of a theo- 
retical model, such as the reliability of the extrapolation 
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FIG. 1: The binding energies of the N = 8, 20, 28, 50, 82, and 126 isotones as functions of the proton number, Z. The results 
obtained from the deformed RMF+BCS calculations with the TMA parameter set (open circle) are compared with the available 
experimental data (solid square) |3S| . 



to the unknown areas. In Fig. 1, we compare the binding 
energies of all the isotones with neutron number N = 8, 
20, 28, 50, 82, and 126 with the available experimental 
data [33. Remarkable agreement between theory and 
experiment can be clearly seen. For the N = 8 isotonic 
chain, the experimental trend is reproduced quite well. 
The only small discrepancy is that the calculated results 
are somewhat systematically larger than the experimen- 
tal data. Nevertheless, it is surprising that a mean field 
model can reproduce such light nuclei so well without fit- 
ting the parameters particularly for this region. For the 
N = 126 isotonic chain, theory agrees well with exper- 
iment around 208 Pb. The results for nuclei with more 
protons begin to deviate from the experimental data, i.e. 
the calculated results are larger than the experimental 
data. The use of other parameter sets, NL3 |3^] and 
NL-Z2 ^(jj does not change this conclusion. 

It is well known that the modern HF or RMF cal- 
culations still can not provide a description of nuclear 
masses comparable to those of most mass tables [4l| . 
This feature is attributed to the limited number of nu- 
clear masses taken into account in the fitting procedure 
of the model parameters. On the other hand, when more 
masses are included into the fitting procedure as what has 
been done by Goriely et al. [12, E3, 03 HH j a description 
of nuclear masses comparable to that of the finite-range 
droplet model (FRD M) l3al has been obtained based on 
the HF and HFB [Hill 111 methods. In Table I, we 
compare the root-mean-square (rms) deviation er of our 
present calculation, with that of the HFBCS-1 mass for- 
mula, and that of the FRDM model for 107 nuclei whose 



experimental (or extrapolated) masses are compiled in 
Ref. H3. 



TABLE I: The mass rms deviation a for 107 nuclei with 
N = 8,20,28,50,82, and 126 from the present calculation 
(RMF+BCS), the HFBCS-1 mass formula (HFBCS-1) 0, 
and the FRDM mass formula (FRDM) jjjjj. 





RMF+BCS 


HFBCS-1 f33] 


FRDM [361 


a 


2.873 


0.957 


0.774 



It is easily seen that the rms deviation of the present 
RMF calculation is still 3 times larger than that of the 
HFBCS-1 mass formula or 4 times larger than that of 
the FRDM model. It should be noted that the rms devi- 
ations for the HFBCS-1 mass formula and for the FRDM 
mass formula shown in Table I are larger than their over- 
all deviations, 0.738 MeV and 0.669 MeV, respectively, 
mainly due to the relatively large discrepancies for those 
light isotones with N = 8, 20, and 28. It is expected that 
a better description of nuclear masses in the RMF frame- 
work be obtained if more masses are taken into account 
in the parameter fitting procedure. 

IV. TWO- AND ONE-PROTON SEPARATION 
ENERGIES 

The two- and one-proton separation energies 

S 2p {N,Z) =B(N,Z)-B(N,Z-2), (4) 



4 




8 10 12 14 

Proton Number Z 



I ' I i I i I ' I 



N=20 



I I I I I I I 



> 20 ; 

0) 15 - 



— ■ — Exp. 

-a- RMF+BCS/TMA 



10 

£ 5 



'0-9 : 

10 12 14 16 18 20 22 24 26 28 30 32 

Proton Number Z 



N=28 



Vv 



- Exp. 
RMF+BCS/TMA 



14 16 18 20 22 24 26 28 30 32 34 36 

Proton Number Z 



50 

40 

> 30 
0) 

5 20 

a 10 



N=50 



H — 1 — h 

o°oo. 



•0°nO 



— ■- Exp. 

b^a q -o- RMF+BcsrrMA 



28 32 36 40 44 48 52 

Proton Number Z 



> 20 

a 15 

•L 10 



V 



N=82 



— ■— Exp. 

-o- RMF+BCS/TMA 



44 48 52 56 60 64 68 72 76 80 

Proton Number Z 



— 2 





N=126 



I ' I i I i I ' I 



- Exp. 
RMF+BCS/TMA 



-----P-----P-- 
b b 



2 84 86 88 90 92 94 96 98 100 

Proton Number Z 



FIG. 2: The two- and one-proton separation energies, S2 P and S p , of the N = 8, 20, 28, 50, 82, and 126 isotones as functions 
of the proton number, Z. The results obtained from the deformed RMF+BCS calculations with the TMA parameter set (open 
circle) are compared with the available experimental data (solid square) |3al . 



S P (N,Z)=B(N,Z)-B(N,Z-1), (5) 

where B(N, Z) is the binding energy of the nucleus with 
neutron number N and proton number Z ', are quite im- 
portant and sensitive quantities to describe nuclei and the 
corresponding shell closures. In Fig. 2, we plot the two- 
and one-proton separation energies for all the N = 8, 20, 
28, 50, 82, and 126 isotones as functions of the proton 
number Z . A larger drop than its neighboring counter- 
parts in these curves is usually interpreted as a (sub)shell 
closure, and the corresponding nucleon number is called 
"(sub-)magic number". In the present work, we mainly 
adopt such a notion to discuss shell closures. 

In Fig. 2, the calculated two- and one-proton separa- 
tion energies are compared with the available experimen- 
tal data [38j . We should mention that some experimental 
data from Ref. [3^| are obtained by "distant connection" 
and that we do not make any distinction between real 
experimental data and extrapolated data. From Fig. 2, 
it is quite clear that the theoretical calculations agree re- 
markably well with the experimental data in the whole 
region from mass number 11 (^Lis) to mass number 218 
(92 Ui2e). 

For the N = 8 isotonic chain, the Z = 8 (lpi/2) shell 
closure is clearly seen from both the two- and the one- 
neutron separation energies, where the shell quantum 
number in the bracket is that of the shell closure. Both 
experimental data and our calculations also show a dis- 
tinct shell closure at Z — 6 (1^3/2)- For the N = 20 iso- 
tonic chain, in addition to the traditional Z — 8 (lpi/2), 
Z = 20 (W3/2), and Z = 28 (I/7/2) shell closures, a 



new (sub)shell closure (1^5/2) is clearly seen at Z = 14 
• 14S120 has been recognized as a doubly magic nucleus 
by a previous shell- model calculation 42], due to a fea- 
ture characteristic of doubly magic nuclei, i.e. it has a 
Op — Oh ground state and its two lowest excited states are 
intruders. For the N = 28 isotonic chain, two neutron 
shell closures at Z = 20 (W3/2) and Z — 28 (I/7/2) are 
reproduced quite well. For the N = 50 isotonic chain, the 
agreement with the available experimental data is very 
good. Although the Z = 28 (I/7/2) and Z = 50 (I59/2) 
shell closures are still unaccessible experimentally up to 
now, they are clearly shown in our calculations. Based 
on the overall good agreement of the N = 50 isotonic 
chain, we believe that our predictions are reliable. 

For the N = 82 isotonic chain, the Z = 50 (I59/2) 
shell closure is reproduced quite well while the Z = 82 
(lft.11/2) shell closure is already in the unbound region. 
We note that there is a new (sub)shell closure at Z = 58 
(157/2)) as shown in Fig. 2. It was argued long time 
ago that there is a whole "plateau" of stability for all the 
even 58 < Z < 70 nuclei in the TV = 82 isotonic chain 
|43l Eij. i.e. the so-called "changing magicities". Our 
calculations do indicate this (sub)shell closure, which is 
also supported by the experimental data, as shown in 
Fig. 2. However, the different deviation trends for nuclei 
with Z > 58 and Z < 58 indicate that some important 
features around Z — 58 in the N = 82 isotonic chain may 
be missed and/or mistreated in the relativistic mean field 
theory. 

For the N = 126 isotonic chain, immediately, we no- 
tice the seemingly big deviations of our calculations from 
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the experimental data. From both the two- and the 
one-proton separation energies, we note that the exper- 
imental Z = 82 (1/111/2) shell closure is underestimated 
somewhat by the relativistic mean field calculations. Our 
calculations also indicate another new (sub)shcll closure 
(1/19/2) at Z — 92 in agreement with Ref. [45j, where 
92 8 Ui26 is studied as a doubly magic nucleus. We no- 
tice that this is a unique phenomenon in the relativistic 
mean field theory while a large-scale shell-model calcula- 
tion did not find this shell closure 0| ■ We also find that 
the use of other often used parameter sets in this region, 
such as NL3 [3^ and NL-Z2 ^cji does not change our 
conclusion. 

The two-proton separation energy becomes negative 
when the nucleus becomes unstable with respect to two- 
proton emission. Hence, the two-proton drip-line nu- 
cleus for the corresponding isotonic chain is the one with 
two less protons than the nucleus at which Si v first be- 
comes negative. In the same way, we can also name the 
one-proton drip-line nucleus. The predicted two- and 
one-proton drip-line nuclei based on our calculations are 



12^85 26^20, 24^20; 

152 V K • 220p.. . 218 
I 1382, oa rUr" 



v , .100 
, 28 iN1 28, 50 



Zn 28 , ^Ni 



Sn 50 



156 
! 74 



w 8 



70 * "82; §4 fui26i 92~Ui 2 6; respectively. It is interesting 
to note that in some cases the two-proton drip-line nu- 
cleus and the one-proton drip-line nucleus are the same 
one, such as i^Mgg in the N = 8 isotonic chain and 
5Q°Sn5o in the N = 50 isotonic chain. In other cases, 
the one-proton drip-line nucleus comes before the two- 
proton drip-line nucleus, which is easy to understand be- 
cause the pairing interaction can increase the stability of 
even-even nuclei compared with its one-neutron less iso- 
tone. Depending on the pairing strength of the specific 
situation, the difference in neutron number between the 
two- and one-proton drip-line nucleus could be two or 
four. More specifically, the difference is 2 in the N = 20, 
28 and 126 isotonic chains, while it is 4 in the N = 82 
isotonic chain. 



V. NUCLEAR RADII 

The difference between neutron and proton root-mean- 
square (rms) radii, 



R. 



(p) 



'n(p)/ 



\l/2 



Jpr 



(P)' 



2 dr 



J Pn( P )dr 



1/2 



(6) 



where is the neutron (proton) density distribution, 
can be very large for nuclei with exotic isospin ratios, thus 
forms the so-called neutron (pro ton) skin [3(]], and even 
neutron (proton) halo [231 124. l2q | . The neutron rms radii 
(empty square) and proton rms radii (empty circle) of the 
N = 8, 20, 28, 50, 82 and 126 isotones are plotted in Fig. 
3, where the experimental proton radii (solid circle with 
error bar) are extracted from the experimental charge 
radii |47| by the formula, 



Since the HFBCS-1 mass formula has achieved an un- 
precedented success in describing nuclear charge radii liSl 
compared with the compilation of Nadjakov et al. |49i |. 
its results are also shown for comparison (open stars). 

In Fig. 3, it is clearly seen that the agreement between 
our calculations and the experimental data is remarkably 
good. For 17CI20 and 22T128! the experimental data are 
somewhat larger than our theoretical predictions. Also, 
the proton radii of two nuclei in the N = 126 isotonic 
chain, §2 8 Pbi26 and |g 9 Bii26, seem to be overestimated a 
bit in our calculations. Some kinks in Fig. 3 are due to 
the sudden changes of deformations of the corresponding 
nuclei (see Fig. 4), including the one at Z = 16 in the 
N = 28 isotonic chain, the one at Z = 25 in the N = 50 
isotonic chain, and the one at Z = 97 in the N = 126 iso- 
tonic chain. It should be noted that the last few nuclei in 
each isotonic chain shown in Fig. 3 are unbound in our 
calculations with separation energies — 5 < S p (2p) < 
MeV (see Fig. 2), thus the sudden increases of the pro- 
ton radii do not necessarily imply appearances of proton 
halos. 

The agreement between our results and those of the 
HFBCS-1 33 mass formula is also very good in general, 
particularly for those nuclei in the N = 50 and 82 isotonic 
chains. The discrepancies for some nuclei with Z w 12 
in the N — 20 isotonic chain are due to the deformation 
effects (see Fig. 4). And the same is true for the nuclei 
with 17 < Z < 30 in the N = 28 isotonic chain, but the 
experimental data seems to be in favor of our results. For 
all the nuclei in the N — 126 isotonic chain, the results 
of the HFBCS-1 mass formula are a bit smaller than our 
results and therefore are closer to the experimental data. 

One more thing we note is that unlike proton radii, 
which go up almost monotonously with increasing proton 
number, neutron radii usually go down first, then go up 
slowly. This phenomenon is more obvious in the N = 8, 
20, 28, and 50 isotonic chains. The underlying reason is 
not difficult to understand. For proton-deficient nuclei 
with small Z, the protons occupy "inner" single-particle 
states with much smaller rms radius than the neutron 
rms radius. Therefore, an extra proton added to these 
"inner" single-particle states attracts the neutrons closer 
to the nuclear center. Consequently, the neutron rms ra- 
dius decreases until Pauli principle prohibits further pro- 
tons from occupying these "inner" single-particle states. 
With more protons added, the opposite becomes true and 
the neutron rms radius begins to increase, though very 
slowly. 



Rt = Rt — 0.64 fni 



(7) 



VI. DEFORMATIONS 

Deformation is another important property of nuclei. 
It also could be used as one of the indicators of magic- 
ity conserving or losing of the corresponding nucleus. In 
our previous work |2(i| , we found that Sn isotopes are de- 
formed in the neutron-rich region, which could be viewed 
as the magicity losing of proton magic number Z = 50. 
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FIG. 3: The rms neutron and proton radii, R n and R p , of the N = 8, 20, 28, 50, 82, and 126 isotones as functions of the proton 
number, Z. The results obtained from the deformed RMF+BCS calculations with the TMA parameter set (open circle) are 
compared with the available experimental data (solid circle with error bar) [i^ and those of the HFBCS-1 mass formula (open 
star) |33|. 



.-■-Exp. j 






-o-RMF+BCS/TMA 






* FRDM 


N= 


8 ] 


-ft- HFBCS-1 






: ■ ■ * 












p, .*/* 






,Q "O-O-^— -&-£> 


-6 


vb-O-0 













2 4 6 


8 10 12 14 16 1 


Proton Number Z 




n=5o ; 




— ■— Exp. 

-o-RMF+BCS/TMA . 




* FRDM 


9' 


-ft- HFBCS-1 


" ■ . « 

^^ss^m^k^ 


- S ftftftft3ftG OGC 

ftft 





24 28 32 36 40 44 48 52 56 

Proton Number Z 



-0.1 
-0.2 - 
-0.3 
-0.4 



N=20 



— ■— Exp. 

-o- RMF+BCS/TMA- 

* FRDM 
-ft- HFBCS-1 



q Q Q 6*.000'*^%" 



8 10 12 14 16 18 20 22 24 26 28 30 32 

Proton Number Z 



N=82 

— ■— Exp. 

-o-- RMF+BCS/TMA 

* FRDM 
-ft- HFBCS-1 
■■>■■ 



52 56 60 64 68 72 76 80 84 

Proton Number Z 



N=28 



- Exp. 

- RMF+BCS/TMA 
FRDM 

- HFBCS-1 



/ O.ft J \ i \ * 

5 V «*, 



20 22 24 26 28 30 32 34 36 38 

Proton Number Z 



N=126 

— ■— Exp. 

-o- RMF+BCS/TMA 

* FRDM 
-ft- HFBCS-1 

bo 



72 74 76 78 80 82 84 86 88 90 92 94 96 98 100 

Proton Number Z 



FIG. 4: The mass quadrupole deformation parameters, 02, of the N = 8, 20, 28, 50, 82, and 126 isotones as functions of the 
proton number, Z. The results obtained from the deformed RMF+BCS calculations with the TMA parameter set (open circle) 
are compared with the predictions of the FRDM model (solid star) |3g, those of the HFBCS-1 mass formula (open star) |33f, 
and the available experimental data (solid square with error bar) extracted from the B(E2 : + — > 2 + ) values |51|. Note that 
the extracted deformation parameter does not mean that the corresponding nucleus is really deformed. 
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However, experimentally, it is difficult to obtain deforma- 
tion knowledge of nuclei directly. One of the most usual 
methods is to extract deformation parameters /3 2 from 
the B(E2 : + — > 2 + ) values by using the assumption 
of deformation [HI]]. Hence, all the nuclei have a finite 
deformation and also the analysis is limited to even-even 
nuclei. Due to this fact, we have to treat the extracted 
deformation parameters with care. For most cases, if the 
extracted deformation parameter is small, the nucleus is 
expected to be spherical. 

In Fig. 4, we plot the mass deformation parameters, 
f3 2 , for all the N = 8, 20, 28, 50, 82, and 126 isotoncs 
as functions of the proton number Z. The experimen- 
tal data are extracted from the B(E2 : + — > 2 + ) val- 
ues given in Ref. The predictions of the FRDM 
mass formula [3(| and those of the HFBCS-1 mass for- 
mula t 33| are also shown for comparison. From Fig. 4, 
it is easily seen that the agreement between the results 
of our calculations and those of the other two methods 
is remarkably good. All these three methods obtain es- 
sentially the same results for all the N = 50, 82 and 126 
isotones. For the N = 8 isotonic chain, FRDM predicts 
igNes, iiNag, and 2 2 Mg 8 to be slightly deformed, while 
in our calculations, these nuclei are spherical. On the 
contrary, the HFBCS-1 mass formula predicts all the nu- 
clei with Z > 10 to be slightly deformed except for 2 |Ss. 
For the N = 20 isotonic chain, both our calculations and 
the FRDM method show that all the nuclei are spheri- 
cal, while the HFBCS-1 mass formula predicts the nuclei 
with Z = 11, 12, 13, and 25 are deformed, which seems 
to be closer to the experimental data. For the N = 28 
isotonic chain, our calculations predict f|S 2 s to be de- 
formed, which is consistent with both the experimental 
data and the predictions of the RHB method |52|, while 
FRDM predicts it to be spherical. These two methods 
also differ in the predictions for the proton drip-line nu- 
clei, iBr 28 , l|Kr 28 , l^ 3 Bk 126 , and || 4 Kr 126 . We should 
note that these nuclei are already unstable against pro- 
ton emission in our calculations (see Fig. 2). On the con- 
trary, the HFBCS-1 method agrees with our calculations 
at both ends of this isotonic chain, but differs a bit in the 
middle: It predicts all the nuclei with 17 < Z < 30 to be 
slightly deformed. That could be the reason why there 
are some discrepancies between our results and those of 
the HFBCS-1 mass formula, as shown in Fig. 3. For 
the N = 50 isotonic chain, the main differences between 
the predictions of these three methods exist in the region 
with 20 < Z < 26. Our calculations show that these 
nuclei are deformed, while both FRDM and HFBCS-1 
predict them to be nearly spherical. 

Except for nuclei with Z < 18 and Z > 30 in the 
N = 28 isotonic chain, and nuclei with 20 < Z < 26 
in the N = 50 isotonic chain, all the other nuclei in our 
calculations are more or less spherical. Although, at first 
sight, the agreement between our calculations (including 
the predictions of the FRDM model and the HFBCS- 
1 mass formula) and the experimental data is not very 
good, we should keep in mind the limitations of the ex- 



tracting method that we obtain the experimental defor- 
mation parameters. Within such a method, even the 
well-known spherical nuclei, g 6 Os and |[JCa 2 o, are pre- 
dicted to be deformed [5l| (see Fig. 4). Therefore, the ex- 
perimental data should be compared with caution. How- 
ever, combined with other experimental knowledge, such 
as the reduced electric transition probability, B(E2), and 
energy, we can decide whether a nucleus is really de- 
formed. The corresponding quantities for JpNes, f|Mg 2 o, 
and gTi 20 are 0.0269(26) e 2 b 2 , 1887.3(2) KeV; 0.039(7) 
e 2 b 2 , 885.5(7) KeV; 0.087(25) e 2 b 2 , 1554.9(8) KeV; re- 
spectively, which indicate these nuclei are really deformed 
|5ll | . The fact that our calculations did not find deformed 
ground-state configurations for these nuclei is probably 
due to the broken rotational symmetry in the relativistic 
mean field theory, which will be discussed below . 

One of the long unsolved problems in the relativistic 
mean field theory is that it cannot obtain the deformed 
ground-state configurations for 32 Mg, which is known 
to be strongly deformed experimentally and other 
N = 20 isotones around Z = 12, including 31 Na and 
30 Ne E3|. On the other hand the AMD method 
[l7| , the shell-model calculations |42t l54| , and the angu- 
lar momentum projected generator coordinate method 
|55| can obtain a deformed ground-state configuration 
for 32 Mg. In the shell- model calculations [I2I IB^j. it is 
demonstrated that the 2p — 2h intruders dominate the 
ground states of 30 Ne, 31 Na, and 32 Mg, therefore lead 
to deformed ground states for these nuclei. While the 
studies of the AMD method and the angular momentum 
projected generator coordinate method clearly show that 
the zero-point energy associated with the restoration of 
the broken rotational symmetry j5|| is indispensable for 
a correct description of the ground state of 32 Mg. 

As mentioned above, we have performed the 
quadrupole constrained calculations for every nucleus in 
our present work in order to obtain its potential energy 
surface and determine the corresponding ground state. In 
Fig. 5, we plot the potential energy surface of 32 Mg as a 
function of the deformation parameter /3 2 . For compari- 
son, the results calculated with the parameter sets NL3 
and NL-Z2 are also shown. For all the three parameter 
sets, an indication of a second minimum around /3 2 ~ 0.5 
is clearly seen (though very "soft" for TMA and NL3), 
which ag rees well with the unprojected AMD an d 
HFB [43] calculations. (Here it should be noted that the 
total energy calculated with NL-Z2 seems to deviate from 
those obtained with TMA and NL3 somewhat, therefore 
does not agree with experiment very well. It could be due 
to the fact that NL-Z2 is aimed at heavy and superheavy 
nuclei but not light nuclei like 32 Mg. Also, we should note 
that TMA and NL3 adopt the same scheme of center- 
of-mass correction, while NL-Z2 uses the so-called mi- 
croscopic center-of-mass correction ^(j •) Therefore, we 
conclude that to obtain the deformed ground-state con- 
figuration for 32 Mg, including 18 Ne and 44 Ti, the broken 
rotational symmetry in the relativistic mean field the- 
ory has to be restored. Since so far, to our knowledge, 



8 




-252 - 



-254 - 

-256 I ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 

-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 

h 

FIG. 5: The potential energy surfaces of 32 Mg as functions of 
the deformation parameter, 02. The solid, dashed, and dotted 
lines correspond to the results calculated with the parameter 
sets TMA, NL3, and NL-Z2, respectively. 

such a calculation is still missing in the RMF method, no 
wonder no RMF calculations have succeeded in obtain- 
ing deformed ground-state configurations for these nuclei. 
We also consider this as one motivation to introduce the 
angular momentum projection method into the present 
RMF model to restore the broken rotational symmetry 
as one of our next works. 

In Ref. [13], the authors argued that their experimen- 
tal high Qp value, 8344±^5 KeV; for 130 Cd is a direct 
signature of an N — 82 shell quenching below 132 Sn. 
While in our calculations, from the point of view of de- 
formations, there is no indications of shell quenching for 
the whole N — 82 isotonic chain, even in the neutron- 
rich side. In our calculations, the ground-state binding 
energy of 130 In, E = 1080.966 MeV (/? 2 = 0.041), to- 
gether with the binding energy of 130 Cd, E = 1072.897 
MeV (/3 2 = 0.0), give a Qp of 8849 KeV, which falls into 
the category of "quenched models" defined in Ref. [l2| . 
Therefore, whether or not the N — 82 shell is quenched in 
the neutron-rich region needs a more careful study both 
experimentally and theoretically. 

VII. CONCLUSIONS 

Nuclei with magic numbers are very important in the 
study of nuclear physics, both theoretically and experi- 
mentally. The localized feature of "magic numbers" has 
received more and more attention due to its importance 
to astrophysical problems and recent experimental devel- 
opments to produce exotic nuclei in the laboratory. In 
the present work, we have adopted the relativistic mean 



field theory to study neutron magic nuclei with six clas- 
sical magic numbers N = 8, 20, 28, 50, 82, and 126. The 
relativistic mean field theory that has proved to be very 
successful in describing many nuclear properties, due to 
its natural spin-orbit description and only a few param- 
eters, is considered to be an ideal model to study not 
only stable nuclei but also exotic nuclei. In order to 
study all those nuclei in a proper way, we must take into 
account pairing correlations and deformation effects si- 
multaneously. For nuclei with odd number of nucleons, 
the blocking effect must also be treated properly [57j . 
All these requirements make the calculation very com- 
plicated and time consuming, so up to now a systematic 
study in the relativistic mean field theory including all 
these neutron magic numbers is still missing. However, 
due to the importance of these nuclei to various physical 
problems, such a study is in urgent need. This work is 
also part of our series of efforts to study the magicity of 
magic numbers in the relativistic mean field theory. 

In our systematic study of the neutron magic isotonic 
chains, we have found a new proton magic number Z = 6 
in the N = 8 isotonic chain. The shell closure feature at 
Z = 14 in the N = 20 isotonic chain is also reproduced 
very well. We have failed, however, to obtain deformed 
ground-state configurations for nuclei around Z = 12 in 
the N = 20 isotonic chain, which have been attributed 
to the broken rotational symmetry. N = 28 has been 
predicted to lose its magicity in the proton-deficient side 
of the isotonic chain, i.e. they are found to be deformed. 
A = 50 seems to conserve its magicity in the proton-rich 
side while lose its magicity in the proton-deficient side. 
This is in agreement with its proton counterpart, Z = 50, 
which is also found to lose its magicity in the neutron-rich 
side. For the A~ = 82 isotonic chain, we have confirmed 
that there is a (sub)shell closure at Z = 58, but the 
different deviation trends in the Z > 58 and Z < 58 re- 
gions seem to imply that some features are missed and / or 
mistreated in the relativistic mean field theory. For the 
A = 126 isotonic chain, both the two- and one-proton 
separation energies are a bit overestimated. Our calcu- 
lations also indicate a new (sub)shell closure at Z = 92. 
The binding energies agree well with experiment around 
208 Pb and become larger than experiment with increasing 
proton number Z. The use of two other parameter sets, 
NL3 and NL-Z2, does not change the above conclusion. 
The predicted two-proton and one-proton drip-line nuclei 
for the N — 8, 20, 28, 50, 82, and 126 isotonic chains are 

20-\/r rr . 46t?„ 44/~V . 587 56t\t: . 100c„ . 156w 

12 iVi g8, 26-^20, 24^ r 20, 30^28, 28 1Nl 28, 50 Sn 50, 74 W 82, 

7 o 2 Yb 82 ; 94°Pui2 6 ; 92 8 Ui2e; respectively. 
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